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We present theoretical description of conduction electrons interacting with a domain wall in fer- 
romagnetic metals. The description takes into account interaction between electrons. Within the 
semiclassical approximation we calculate the spin and charge distributions, particularly their modi- 
fication by the domain wall. In the same approximation we calculate local transport characteristics, 
including relaxation times and charge and spin conductivities. It is shown that these parameters are 
significantly modified near the wall and this modification depends on electron-electron interaction. 
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I. INTRODUCTION 

It is generally believed that domain walls modify sig- 
nificantly all the transport properties of ferromagneLin 
metals. Early theoretical analysis of this phenomenoro'El 
were stimulated by magnetotransport measurements on 
single-crystal Fe-wiskers. Recently progress in controlling 
magnetic state of nanostructures enabled observation of 
a direct correlation between domain structure and trans- 
port properties.u For example it is possible to extract 
the contribution to resistivity due to a single domain 
wall. In addition, the discovery of Giant Magnetore- 
sistance (GMR) in magnetic multilayers, which is con- 
nected with reorientation of the magnetic moments of 
neighboring magnetic layers from antiparallel to parallel 
alignments, renewed the interest in domain wall resis- 
tivity. To some extent, the domain wall plays a similar 
role as the nonmagnetic layer separating two ferromag- 
netic films in a sandwich structure or in a multilayer, and 
therefore can be expected to lead to magnetoresistance 
effects similar to GMR. Indeed, there is a growing exper- 
imental evidence of a large magnetoresistancCj-due to a 
domain wall in ferromagnetic nanostructures.EnJ This, in 
turn, led to growing interest in theoretical understand- 
ing of the behavioE, oL electrons coupled to a ferromag- 
netic domain wall.Enlil Moreover, progress in nanotech- 
nology made also possible to study electric current flow- 
ing through a narrow contact between two ferromagnetic 
metals (point contact), where a constrained domain wall 
is createdj-in the antiparallel configuration. Such a do- 
main wallcJ has a significant influepce on the transport 
characteristics of the point contact Jil 

It has been shown experimentallyllJ that the presence 
of a domain wall can either increase or decrease elec- 



trical resistance of a system. This intriguing observation 
stimulated theoretical works on understanding of the role 
of a domain wall in transport properties. Levy et alE2l 
developed a semiclassical model based on the mixing of 
spin-majority and spin-minority transport channels by 
the domain wall. This mixing results in increase of the 
electric resistance due to, the presence of a wall. On the 
other hand, Tatara et alB found a negative contribution, 
which is due to destruction of the weak localization cor- 
rections to conductivity by the domain wall. Another 
model which may lead either to positive or negative con- 
tribuiion of a wall to resistivity was developed by Gorkom 
et al£3 The key point of this model is the fact, that the 
wall can lead to redistribution of the charge carriers be- 
tween spin-majority and spin-minority channels. The do- 
main wall contribution to resistivity depends then on the 
ratio of spin-majority and spin-minority relaxation times. 

In this paper, we consider electrons in a ferromagnetic 
metal, which interact with a domain wall. The descrip- 
tion includes interaction between electrons, and therefore 
we use a self-consistent analysis to describe charge and 
spin distributions, as well as their modification by the 
domain wall. Using the Green's function technique, we 
calculate the electron relaxation times in the quasiclassi- 
cal approximation. Apart from this, we also calculate the 
local charge and spin conductivities. These transport pa- 
rameters are shown to be significantly modified near the 
domain wall, which may give rise to new effects. 

The paper is organized as follows. In Section 2 we de- 
scribe the model as well as the transformation used to 
replace the system with inhomogeneous magnetization 
by a system magnetized homogeneously. In Section 3 we 
present the transformed Hamiltonian, generalized by in- 
cluding selfconsistent fields related to electrostatic and 
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magnetic interactions. A semiclassical solution of the 
resulting Schrodinger equation for electrons is also pre- 
sented there. Scattering from the wall, in the Born ap- 
proximation, is calculated in Section 4. In Section 5 we 
calculate, within the quasiclassical approximation, both 
the spin and charge distributions in the vicinity of the 
domain wall, as well as the corresponding contributions 
generated by the wall. Local relaxation times are cal- 
culated in Section 6, whereas the local charge and spin 
conductivities are calculated respectively in Sections 7 
and 8. Final conclusions are provided in Section 9. 



II. MODEL 



Consider a general case of a ferromagnet with a nonuni- 
form magnetization M(r). The one-particle Hamiltonian 
describing conduction electrons locally exchange-coupled 
to the magnetization M(r) takes the form 

1 9^ 

Ho^-—4'lj-^iJc.-J^io-^p-M{r)i;^, (1) 
2m OT^ 

where J is the exchange parameter, ipa and ipl^ are the 
spinor field operators of electrons, a = {ox,<Jy,<7z) are 
the Pauli matrices, and we use the units with h = \. 

The model Hamiltonian (1) will be used to describe 
electrons interacting with a domain wall in a ferromag- 
netic metals or in semiconductors. The domain wall will 
be modeled by a magnetization profile M(r). For the 
sake of simplicity we shall assume that |M(r)| = const. 
We can then write 



JM(r) = Afn(r), 



(2) 



where n(r) is a unit vector field to be specified later, and 
M measured in the energy units includes the parameter 
J. 

The first step of (war analysis is to perform a local uni- 
tary transformationB 



^^T{v)^, Tt(r)r(r) = l, 



(3) 



which removes the nonhomogeneity of M(r), that is T(r) 
transforms the second term in Eq. (1) as 



■0^ (T • n(r) ip ip^UzTp- 



(4) 



The transformation matrix T(r) must then obey the con- 
dition 

Tt(r)<T.n(r)T(r) = (T,. (5) 
Explicit form for such a T(r) is given byEl 

„, N 1 / I i-T (r) — n,. (r) cr,, \ , , 

T{v) = ^[^l + nz{v)+i . (6) 

V2 \ ^l + nz{v) 



The transformation (3), (6) can be applied not only to 
a simple domain wall, but also to other types of topo- 
logical excitations in ferromagnetic systems, for instance 
helicoidal waves, skyrmions, and others. 

Applying the transformation (6) to the kinetic part of 
the Hamiltonian (1) one obtains 

^'£V'-V'^(|: + A(r))%, (7) 

where the non-Abelian gauge field A (r) is given by 

A(r)=Tt(r)|-T(r). (8) 

According to Eq. (6), the gauge field A(r) is a matrix in 
the spin space. 

Let us consider now a more specific case of a domain 
wall in a bulk system. Assume, the wall is transla- 
tionally invariant in the x-y plane: M(r) M(z) and 
n(r) n(z). For a simple domain wall with M(z) in the 
plane normal to the wall, one can parametrize the vector 
n(z) as 



n{z) — {si\np{z), 0, cos {p{z) ) , 



(9) 



where the phase ^p{z) determines the type of a domain 
wall. The transformation (6) is then reduced to 

V2 y Vl + cosip{z) J 



and the gauge field assumes a simple form 



A(z)= ( 0, 0, --a, ^'{z) 



(10) 



(11) 



where (p'{z) = dip{z) / dz. 



V-f 



FIG. 1. Schematic picture of the magnetization orienta- 
tions near the domain wall. 
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Taking into account Eqs. (7) and (11), one can write 
the transformed Hamiltonian in the form: 



Ho 



1 mf3^{z) 



-la,, 



2 

iay I3{z) 



d_ 

dz ' 



where 



P{z) 



2m 



(12) 



(13) 



For a slowly varying smooth function ip{z) (thick domain 
wall centered at z=0 with width L), the perturbation 
due to the domain wall is weak, and close to the center 
of the wall, \z\ <C L, the parameter (3{z) can be treated 
as a constant. Such a model domain wall with a constant 
parameter (3 was analyzed in Ref . [ ^| . 

The description given above is quite general and may 
be used for various models of the domain wall. If we 
assume the domain wall in the form of a kink (Fig. 1), 
then 



tanh {z/L), 



and the parameter /3(z) is given by 



AmL cosh {z/L) 



(14) 



(15) 



III. SEMICLASSICAL APPROXIMATION 

In bulk magnetic metals like Fe, Ni or Co, the width 
i of a magnetic domain wall is usually much larger then 
the electron Fermi wavelength A^. In such a case applk 
cation of a semiclassical approximation is well justifiedJll 
The dominant perturbation from the domain wall is then 
described by the term with /3 (d/dz) in Eq. (12), since it 
is of order oi pkp- The term proportional to 0^ is smaller, 
while the term including /3'(z) is of the order of (3/L and 
therefore can be neglected. 

For the sake of self-consistency, we will include now 
the Coulomb interaction of electrons, which allows cor- 
rect description of charge accumulated at the wall. The 
point is that the wall can give rise to some excess charge 
locally breaking electrical neutrality, as will be described 
in more details later. This rsffect was not taken into 
account in previous analysis.^ On the other hand, the 
renormalization of the chepical potential forbidding the 
formation of excess chargal3 may be an over-estimation 
of the Coulomb repulsion. 

The Coulomb interaction will be taken into account 
via the coupling term 

H,nt = ljd''r [^/.t (r, t) ^(r, t) - m] F(r - r') 

X [^t(i.',t)^(r',i)™no] , (16) 



where y(r — r') = e^/ |r — r'| is the bare Coulomb inter- 
action and riQ is the mean electron density in the bulk. 
Using an auxiliary scalar field '/'(z) we can incorporate 
the interaction by adding to the Hamiltonian the follow- 
ing term (Appendix A): 



where the field 
equation 



(z) Ip'' -0 



(17) 



(z) is determined by the saddle-point 



dz2 



(18) 



with ( ) denoting the ground state average. This makes 

the solution self-consistent, and the field (j){z) is the 
mean-field electrostatic potential in the presence of the 
wall. The use of the differential saddle-point equation for 
(f>{z), Eq. (18), makes the problem more complicated due 
to the nonlocality, but allows to describe correctly the 
screening effects associated with a spacial distribution of 
charges in the vicinity of the domain wall 

To include the spin-dependent interaction, we will use 
a simpler formalism. More specifically, we introduce the 
contact interaction term in the form 



So 



(19) 



where gs is the corresponding coupling constant. Choos- 
ing So as the spin density far from the wall guarantees 
that this interaction vanishes when there is no domain 
wall. This means that the effects due to magnetization 
of the conduction electrons in a system without domain 
wall are included by the parameter M in the one-particle 
Hamiltonian. The effect of a domain wall is then to mod- 
ify the internal magnetization, resulting from a redistri- 
bution of the spin density. The effects due to interaction 
(19) can be taken into account by adding to the Hamil- 
tonian a new term. 



(fir mz{z) Gz "0, 



(20) 



where the internal magnetization field mz(z) is deter- 
mined (viz. Appendix A) by the saddle-point like equa- 
tion, and is of the form 



Tnziz) = ~gl ((^^2-0) - So) 



(21) 



Thus, the Schrodinger equation for an electron of en- 
ergy e in the fields 0(z) and m^(z) reads: 



2m dr^ 



mP'^{z) 



+i<Ty (3{z) 



dz 



[M - mz{z)\az 



z) - e -0 = 0, 



(22) 
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where the fields 0(z) and m^Jz) have to be determined 
self-consistently via Eqs. (|lqj2l| ). Equation (22) has the 
following semiclassical solutions (i = 1,2) 



exp(±i q • p) 



[sUz)+P^{z)kf{z)]^/ny'{z) 
zki I ki{z) dz 



^iBkAz) 



(23) 



where p = {x,y), q is the momentum in the plane of 
the wall, the wavevector components normal to the wall 
(along the axis z) are given by 



±2m [M^{z)+P^{z)k^{z)]^^^ , 



' " 2to + 2 2m ■ 

and k{z) and Mr{z) defined as 

K^{z) = 2m [e + (t){z 



Mr{z) = M - m^{z). 



(24) 



(25) 



(26) 



(27) 



For clarity of notation we omitted the z label of ki,2{z) 

[kz-iA^) = ^1,2(2:)]- 

There is no reflection from the wall in the semiclassi- 
cal approximation. It should be noted, however, that for 
electrons moving nearly parallel to the wall (with very 
small z-component of the momentum), there is a reflec- 
tion from the wall since for such electrons the classical 
motion through the barrier is impossible. We do not con- 
sider here this effect since its contribution is very small. 



IV. SCATTERING FROM THE WALL (BORN 
APPROXIMATION) 

For the case of not too thin domain walls, the term 
proportional to I3{z) d /dz can be treated as a small per- 
turbation and scattering from the wall can be evaluated 
within the Born approximation. The matrix elements of 
the ikz t) —> {k'2 I) spin-flip scattering in the plane wave 
basis is given by 



Vk^y 



-/3(z)— le 
dz 



ik. z 



dz. 



Using (15), we obtain 

inkz f°° exp[—i(k'z — kz)z]dz 



cosh {z/L) 



(28) 



(29) 



After integrating by parts, this integral can be presented 
as 



_ 1 f°° e~*P^ dx 
Re lim e — - — plm / — 

x^oo 2 Jq 1 -|- e ^ 



where p = {k'^ — kz)L/2. The integral in the last term 
can be calculated using the series expansion of the de- 
nominator in the integrand, and we obtain: 



Vk^k' 



inkz 



m 



1 



i-iy 



n— 1 



Now we use the known representatioruHl 
csch(z) = - + 2z V • ' 



and we find finally 

in^kzik'^ — kz)L 



Vk^k'^ = 



Am 



csch 



(30) 



Correspondingly, the probability of backscattering (fc^ 
-kz) is 



Wback = 2n\Vk^^-kA = 



2to2 



csch^(7rfczL). 



For kzL 3> 1, from the last equation we find 
2^5 k^L^ 



Wba 



ck 



-2-Kk^L 



(31) 



Thus, the probability of the backscattering with simul- 
taneous Sisin-flip vanishes exponentially in the limit of 
kpL ^ l.EJ Spin-conserving backscattering is determined 
by the term proportional to 0^ in the Hamiltonian (22). 
In the first approximation this term can be neglected as 
it is smaller than the term proportional to I3{z) d/dz. 

The question arises, whether the Born approximation 
gives correct results in the problem under consideratinn. 
There are two general conditions for its applicabilitycj 



\U{z)\<. 



lL2 



or \U{z)\ < 



mL ' 



(32) 



where U{z) is the scattering potential. In the first case 
the Born approximation is good for arbitrary electron en- 
ergy, whereas in the second one it is good only for fast 
electrons. Therefore, if we choose the limit kpL ^ 1 then 
|J7(z)| ~ (3kz ~ (kz/mL), and none of the conditions is 
satisfied. In the opposite case of a small domain-wall 
width, kpL ^ 1, we have |J7(z)| ^ 1/mL? and the Born 
approximation is not justified again. Thus, the Born ap- 
proximation can be used only for rough estimations. In 
the case under consideration, kpL 3> 1, it shows that the 
usual scattering from the wall is exponentially weak. 
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V. DISTRIBUTION OF SPIN AND CHARGE 
DENSITIES (SEMICLASSICAL APPROACH) 

In the framework of the semiclassical approximation, 
one can calculate the local charge and spin densities in 
the vicinity of the wall, as well as the distribution of 
charge and spin currents. As follows from Section III, 
this can also be done taking into account the Coulomb 
interaction self-consistcntly. 

The equation for the Green function with the term P{z) 
weakly dependent on z, 



e + 



2m dr'^ 



+ (i>{z)+Mr{z)a, 



d 



-lay /3(z)^ + fij G, (r, r') = 6 {v - v') , 

has a quasiclassical solution (fc^i <C 1), 

e - Ek - Mr{z) (Jz - kz (i{z) ay + 



(33) 



{e - £kT + fJ,r + iS sgne) (e - Skj. + fJ,r + iS sgne) 

(34) 

where the following notation has been used: 

q' + k'z 



£k = 



2m 



£kT,i=^kT[M2(^) + fc,^/?'(-2)] 



ij,r = ^J- ^ ^(P{Z), 



,1/2 



(35) 
(36) 

(37) 



and n is the chemical potential. Equation (36) de- 
scribes the energy spectrum in the spin-up and spin-down 

branches, where for the sake of notational simplicity we 
droped the z-dependence of £kt,i and jJ-r- In what follows 
we also drop the ^-dependence of Mj. and f3. 



Note, that the square root in Eq. (36) contains contri- 
butions due to spin-mixing caused by the perturbation 
(Tyfi [d / dz). Hence, what we call spin-up and spin-down 
branches of the spectrum (labeled as t and [ in Eq. (36)) 
are in fact the eigenvalues corresponding to the wavefunc- 
tions with mixed up and down states. Correspondingly, 
the Green function (34) has poles at both Fermi surfaces 
with k = kpi and k = kpi in diagonal and non-diagonal 
components. 

Using the Green function (34), one can calculate the 
spin density distribution in the presence of the wall. The 
real spin density distribution, i.e., transformed back to 
the original basis is given by the formula 

/dp (V^lf 
2^j2nr^^^'^''^^'^^'^^'^- ^^^^ 

To obtain this expression one should use the inverse of 
the transformation (3), i. e. G^TGT^. Using Eq. (5), 
one can rewrite Eq. (38) as 



„(.).s(.) = -iTry-^a.G.(k). 



(39) 



Since the projection of s on the plane perpendicular to 
n vanishes, we can write the spin density as 



s(.) = -^„(.)Tv/||^..G.(k). 



(40) 



Substituting the Green function (34) into Eq. (40), we 
find 



s(^) 



n(.)(/ -/ 



d^k 



1 



(27r)3 y/TTfcpVW 

(41) 



After evaluating these integrals, we find the spin density 
distribution 



J 



s(^) = 



Mr n{z) 



-'^{M^ + k%,P^f'+'-^{M^ + k%,P^f' 



2/3 



-I- \2mnr + I I arcsmh -j^ arcsmh ^ ) -|- 2mp {kpt + kpi) 



(42) 



r 



where 

kp^^l = 2mnr + 2m^/3^ ± 2m (2m/Xr/3^ + m^P^ 



M, 



(43) 

The accumulation of the spin density at the domain 
wall is 



where So is the spin density in the limit of /3 = 0. For 
2Nji/2 small /3, i.e., for very smooth magnetic wall and up to 
second order of (3 this reads: 



As(2;) = 3(2;) - So, 



. , , 11(2) r 2 /7 . X 
(44) ^s(2;) = [m [kp^ - kpi) 
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7.5 _ 7,5 

10M2 



(45) 



The sign of the factor in the sqnare brackets of Eq. (45) 
depends on material parameters. 

Charge density distribution can be calculated in a sim- 
ilar way, 



After calculating the integral (46) we find 



(46) 



Wc can linearize Eq. (18) in (t>{z) assuming that the 
domain wall does not change significantly the electron 
density, i. e. for ^ ^ |0(2^)|- Hence, after expanding p{z) 
in (j){z) and Fourier transforming over z, we can write 
Eq. (18) as 



(^2 + kI) <j){q,) = -47re2 Ap(g^) 



(48) 



1 /2 

where kq = (47re2 vq) is the inverse screening length, 
vq = dp/dn is the thermodynamic density of states, and 
Ap{qz) is the Fourier transform of 



1 

4^ 



2mp,r [kp] + kpi) 



1.3 I u'A 



mM^ . , kFi/3 mM?: . , kpil^' 
H — arcsmh^r4 — arcsmh- 



Mr 



Mr 



(47) 



Now we can use this expression in Eq. (18) to determine 
electrostatic potential (j){z) and (ip'^'tp) = p{z). 



z)=0 ■ 



(49) 



Using Eqs. (27), (36), (37) and (43), wc find that for 
P ^ the accumulation of charge, Ap{z) = p — p{4i = 0) 
is 

(50) 



This value of Ap{z) is the accumulated charge in the ab- 
sence of Coulomb repulsion. 
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FIG. 2. Distribution of Ap(2:) near the domain wall for 
different values of the spin coupling constant gs- This de- 
pendence corresponds to the charge accumulation in the ab- 
sence of Coulomb repulsion. The distribution of charge in the 

presence of Coulomb interactions, Ap{z), is related to Aj5{z) 
by Eq. (52). The coupling constant gs is given in units of 
ergV^cm^/^. 



30 



25 



o 15 

X 



<l 



10 









0.5x10 




1.0x10"'^ 

1.5x10 

2.0x10"" 


/ \ 
/ \ 
/ . \ 
; \ 
/ ■' ,N '■. \ 




/, /'^^^■,\ 

I ■// \\- \ 

\v\ 

'7/ 

/.-,? \\ 
r'l \\ 




li \ 
i-l \ 
i-i \ 
/■/ ^ 
/■/ 
// 

// 


\ 

\\ 

\ 
\ 


// 
// 


\ 


-2 -1 


1 2 


z/L 





FIG. 3. Distribution of the excess spin density near the 
domain wall for different values of the spin coupling constant 
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Due to (18), the real distribution of charge Ap{z) is 
related to Ap{z) by the relation which in Fourier space 
has the form: 



Tz + i^o 



(51) 



Then, if the characteristic lengh of the domain wall is 
large, kqL ^ 1, we obtain 



^ , ^ 1 cPAp{z) 



(52) 



and, finally, using Eq. (47), we find the distribution of 
charge 



Ap{z) 



m 



d(3{z) 
dz 



(53) 



The results arc presented in Figs. 2 and 3 for differ- 
ent values of the spin coupling constant 9,,. In the cal- 
culations we take the Fermi energies ep'j = 3 eV and 
spi = 2.5 eV in the bulk, and m = 4 mo, where mo is 
the free electron mass. The figures demonstrate how the 
spin coupling constant gs affects both the spin accumu- 
lation (Fig. 3) and Ap{z) (Fig. 2). This effect is a result 
of self-consistency, because by controlling the magnetic 
density one modifies the magnetic wall, and this in turn 
influences the electron density. 

In view of Eq. (49), the function Ap{z), presented in 
Fig. 2, is not the excess charge accumulated at the wall 
but an auxiliary function corresponding to the condition 
(j){z) = 0. The real distribution of accumulated charge, 
Eq. (52), is presented in Fig. 4 for different values of the 
coupling constant gs- This figure demonstrates that the 
integral of Ap(z) over z is zero due to the conservation 
of electric charge. The characteristic length of the charge 
distribution is determined by the characteristic thickness 
of the domain wall. 



We performed numerical calculations of the charge 
Ap{z) and excess spin density Asi{z) = As{z) ■ n(z), 
using the set of equations (42), (47), (49), (27), (37), 
(43), and (21). 
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FIG. 4. Distribution of the charge density near the domain 

wall for interacting electrons, calculated for different values 
of the spin coupling constant Qs- Due to the screening effects, 
the integrated charge accumulated at the wall vanishes. 



VI. IMPURITY SELF-ENERGY 

In this Section we shall take into account the scatter- 
ing of electrons from impurities. The simplest choice are 
impurities with a short-range scattering potential, which 
scatter similarly both spin-up and spin-down electrons. 
Let the matrix element of the scattering potential of a 
defect be Vq. The self-energy operator in the Born ap- 
proximation is 

ne) = VS /-— ^G,(k). (54) 



(2^) 

After integrating the Green function given by Eq. (34), 
we get 



S(£) 



I sgne 



diag(l/T|, 1/t|) 



(55) 



where the momentum relaxation times for spin-up and 
spin-down electrons are 



27r 



Mr . .kF]P Mr . ,k 

H — TT arcsmh — — — arcsmh — , , 

13 Mr P M, 



(56) 



niz) 



m Vi^ 
2tt 



{kpi + kpi 



Mr 

'T 



arcsinh 



kpT P M, 



Mr 



+ 



P 



arcsinh 



kpi P 

Mr 



(57) 



The difference in scattering times is due to a difference 
in the density of states at the Fermi level for spin-up 
and spin-down electrons. The formulae (56) and (57) 
take into account variation of scattering times near the 
domain wall. These correlations have been neglected in 
previous works. 
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VII. LOCAL CONDUCTIVITY 

The general formula for local conductivity (without lo- 
calization corrections), when an electrical field is applied 
along the axis z, has the following form 

g2 r ^3j, 

x{k,-ml3ay)G^, (58) 

where the gauge potential A(z), defined by Eq. (11), is 
taken into account, and the retarded (R) and advanced 
(A) Green functions are both evaluated at the Fermi 
level, 



G«'^^G_f!l{l(k) 

_ -ck - MrUz - kzPUy + /ir ._„s 

(-EkT + /^r ± V2tt) (-£ki + Mr ± V2n) ' 

Owing to the terms containing /3 in the Green func- 
tions, and to the ^-dependence of £fe-f,j,, jJ-r and M^, the 
local conductivity a^z is a smoothly varying function of 
z. Given the azz{z)^ the resistivity of a sample with a 
domain wall can be found by (adding resistivities), 

R-j^y (60) 

Using Eqs. (36) and (37), we find the local conductivity 
in the form 



<7zz = 



27r^m 



n ( -f^ + ™^/^^ ^^-T - rn^Mr/3 arctan 



+ T| [ ^ + kpi- m^Mrl3 arctan 



■ (61) 



One should note that the dependence on j3 enters here not only explicitly, but also through the parameters Tf, t^, 
kpi and kpi- 

The description of the domain wall in terms of local conductivity is justified when L ^ U where I is the electron 
mean free path. For such a smooth domain wall, there is no electron scattering from the wall but the local conductivity 
is changed. The system with a domain wall is macroscopically inhomogeneous, and thus the electric field in the vicinity 
of the domain wall is inhomogeneous when a bias voltage is applied. 



VIII. SPIN CURRENTS AND LOCAL SPIN CONDUCTIVITY 



The spin-current density in the untransformed basis has the form derived in Appendix B: 

de (i^k 



(k). 



(62) 



Suppose the spin current is induced by an electromagnetic field with vector potential Aem acting on both up and 
down spin components. Then, using Eqs. (11), (13) and (62), we obtain for the up and down spin conductivity 



-j^^ [fe - rnfiay ± (n^cr^ - nxOx)i}iz - mf3ay] {kz - mf3ay) Gy. . 



47rTO2 

The result of calculation can be presented in the form 



T I 1 C 

2e 47r^m 



n 



(m/3^ =F Mrnz) {rnkpi 



+{2m0' T MrHz) [j^^M^ + 4^/32 - ^ arcsinh 



— - — arctan — — 

(3 Mr 

_arcsmh^ 

" ^2 „r ^ f , mMr kp-sfi 
[mfi =F MrUz) I m/cFT a — arctan 

\ fj IVlj' 

-(2m/3^ T MrUz) + k%,(3- - ^ arcsinh ^ 



(63) 



(64) 
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In the limit of /? ^ we obtain 

<^lz = TT^zz ± COS z — —V - -^T^ • 65 

The spin conductivity (64) and (65) describes a response 
in the form of up and down spin currents to the electric 
field, associated with the vector potential Aem • 



IX. SUMMARY AND CONCLUSIONS 

We have described behavior of conduction electrons 
interacting with a magnetic domain wall in ferromag- 
netic metals. In the description we used a realistic model 
which includes the Coulomb interaction and screening 
effects. Within the semiclassical approximation we cal- 
culated self-consistently the equilibrium charge and spin 
distribution in the presence of a domain wall. We showed 
that this distribution is significantly modified by the wall. 
We have also calculated the local transport characteris- 
tics, like relaxation times and charge and spin conductiv- 
ities. 

Our approach applies to the linear response regime, 
and therefore such nonequilibrium phenomena like spin 
accumulation at the wall due to flowing current axn not 
taken into account. In a recent paper Ebels et al B ob- 
served large magnetoresistance due to a domain wall and 
attributed it tp-jthe spin accumulation. On the other 
hand, SimanekE3 showed that spin accumulation is par- 
tially suppressed by spin tracking and cannot explain 
such a large magnetoresistance. 

The approach used in Ref. [ |22| is based on the ki- 
netic equation for the Wigner function and takes into 
account nonlinear effects, particularly those due to spin 
accumulation. Such effects were not included in our de- 
scription, since we analysed linear response regime only, 
which is determined by equilibrium characteristics. How- 
ever, we took into account interaction between electrons 
and showed that this interaction can singificanly modify 
the influence of the magnetic wall on transport proper- 
ties. The local transport characteristics were described 
by few parameters characterizing the domain wall. Vari- 
ation of the local conductivities at the wall may lead 
to several new effects. For instance, one may expect the 
Peltier effect at the domain wall. To our knowledge, such 
an effect has not been studied yet. Some other interest- 
ing phenomena may be related to the spin dependent 
coupling, described by the parameter gs ■ 
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APPENDIX A: CHARGE AND 
SPIN-DEPENDENT INTERACTIONS 

To include into the one-particle Hamiltonian the cor- 
rections due to the electron-electron interactions, we use 
a self-consistent mean-field approximation. We can write 
down the partition function Z of our system in form of 
functional integral over spinor fields 

Z = J Dil}'^{Y,t) DiP{r,t) exp (^i J d^rdtL^ , {A.l) 

where Lg is the Lagrangian density: 

Lo = ^t(r,i) (^,^_i7^^(r,t). (A2) 

The contribution to the Lagrangian from the term Hint 
is 

Unt = -\ Spr{t) Vrr' S Pr' (t) , (A.3) 

where T4r' is an infinite matrix with the elements V{r—r') 
and Spr is a vector with elements Sp{r) = ip^ir, t) ip{r, t) — 
riQ. Since the Coulomb interaction is instantaneous, both 
Sp(r, t) are taken at the same time t. 

We use the Hubbard-Stratonovich method enabling de- 
coupling of the interaction term. It gives us the addi- 
tional integration over field 4'{v) in the partition function 

Z = j DtljHr,t) Dij{v,t) D(j){v) exp j cPrdtL 



where 



(^.4) 



L = Lo + 0r <5/0r + ^ (/-r V^^} 0r' . (A.5) 

and, 0r is a vector constructed of the elements 4>{v). 
The mean field approximation corresponds to the saddle- 
point solution for 0(r) 



5L 

5<t>{Y) 



= 0, 



which gives us 



0r' + 5p, = 0. 



(A.6) 



{A.7) 



After Fourier transforming and using ^(q) — 47re'^/g'^, 
we obtain 



<^(q) + 4^6^ 5p{q) = 0, 



(A.8) 
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and, coming back to the r - space, we get the Poisson 
equation for the scalar potential 



+Tt(r)a,T(r) 



d_ 
dr 



■A(r) 



{B.5) 



(A.9) 



In the case of point-like interaction with V{v — t') = 
g 6{r — r'), using the same formalism we obtain the 
saddle-point equation in the form 



(A(r) + g Sp{r) = 0. 



(A.10) 



APPENDIX B: SPIN CURRENT DENSITY 

To find the expression for the spin current, we add to 
the Hamiltonian of Eq. (1) an auxiliary vector potential 
Aem{t) acting only on the spin- up states. It produces in 
the kinetic part of the Hamiltonian the following term 
(in the untransformed basis) 



Hkin = -■^V'^r,*) "5 ^ Am 

2m \or c 2 



After expanding over Aem {t) , we find the linear in Aem {t) 
term in the Lagrangian density 

AL = -H,,r. = -^^r^f^ V^(r,t) (1 |- V(r,t). 

2mc OT 

{B.2) 

The transformation (3) changes it to 



AL = 



+Tt(r)<7,T(r)(-+A(r] 



The corresponding operator of the spin-current density 
can be found by variation 



JT = 



c 6L 



e 5Aem{t) ' 



(BA) 



which gives us finally 
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